Thermal Conductivity as a Probe of Quasi-Particles in the Cuprates. 
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S. A. 



In underdoped YBa2Cu30 x (x = 6.63), the low-T thermal 
conductivity n xx varies steeply with field B at small B, and 
saturates to a nearly field- independent value at high fields. 
The simple expression [1 +p(T)|B|]~ 1 provides an excellent 
fit to Kxx(B) over a wide range of fields. From the fit, we ex- 
tract the zero-field mean-free-path, and the low temperature 
behavior of the QP current The procedure also allows the 
QP Hall angle 9qp to be obtained. We find that 9qp falls 
on the 1 /T 2 curve extrapolated from the electrical Hall angle 
above T c . Moreover, it shares the same T dependence as the 
field scale p(T) extracted from k xx . We discuss implications 
of these results. 

1. Introduction 

Thermal conductivity is potentially a very useful probe 
of the quasi- particle excitations in the superconducting 
state of the cuprates because it is capable of detecting 
the quasi-particle (QP) current in the bulk |l|-Q . In ad- 
dition, measurements of its field dependence may yield 
quantitative information on the QP mean free path. At 
present, this seems to be the best way to investigate 
the low-lying excitations of the condensate. However, 
the task of disentangling the QP current from the larger 
phonon current in the cuprates poses a difficult problem 
for experiment. 

We report recent experiments in which the direct sepa- 
ration of the QP current is achieved by detailed analysis 
of the field dependence of the longitudinal conductivity 
k xx . This line of approach was motivated by the obser- 
vation of plateau features in high-purity P^S^CaC^Os 
(Bi 2212) §]. The existence of the plateaus at low T 
(where dK xx /dH = 0) implies that, in the cuprates, vor- 
tices are essentially transparent to the phonons. Exten- 
sions of these measurements to underdoped YBa2Cu30 x 
(YBCO) reveal that this result may be a rather general 
feature of the cuprates in the clean-limit. This seems to 
us a significant finding since it allows a direct separation 
of the QP current by the application of an intense field. 
In addition, we find that the zero-field mean-free-path £o 
of the QP may be estimated to within a factor equal to 
the vortex scattering cross-section <j tr . 

The isolation of the QP current allows more specific in- 
formation to be extracted from the thermal Hall conduc- 
tivity K xy (Righi-Leduc effect) With the electronic 
current independently determined, we may now obtain 
the Hall angle tanfgp. The QP Hall angle uncovers a 
number of interesting features which we discuss below. 

2. Experiment 



Measurements of the mixed state of the cuprates 

have been reported by several groups. |}|-^) Experiments 
in intense magnetic fields H are complicated by prob- 
lems such as cleaving of the crystal by the large torques 
generated. The most serious problem, however, seems to 
stem from the field sensitivity of the thermometers. At 
the resolution needed, the field dependence of the sensors 
is serious (in thermocouples, moreover, the field sensitiv- 
ity is also history dependent). Previously, we employed a 
bridge-balance method to get around the field- sensitivity 
problem ||. 



14 



12 



10 



E 



o 



AAA A A A A A ^ 0.0 



V 

> 

ft 



M 



A A A A A A A AAA 



-0.2 



D) 

-0.4 E 
CD 



-0.6 







20 



100 



40 60 80 
T(K) 

FIG. 1. The T dependence of the in-plane thermal con- 
ductivity k xx (Q) in zero field (open symbols) in underdoped 
YBCO with x = 6.63 and T c = 63 K. The solid symbols are 
the phonon conductivity Kb estimated by the fit of the ther- 
mal magnetoconductance to Eq. |l| The dc magnetization of 
the crystal is shown as open triangles (H — 5 Oe). 

In our present approach, we adopt a single- heater, two- 
sensor method in which the temperature difference 8T 
between the ends of the sample is detected by two closely 
matched resistive sensors (cernox). The thermal gradient 
— VT is applied in the ab plane, and H is parallel to c. 
The temperature is regulated (with a third base-cernox) , 
and measurements are taken after waiting about 10 min. 
for the field to stabilize. The readings of the two sen- 
sors are recorded at three values of the heater current 
1 = 0, 0.4, 0.6 mA (typically). The 1 = readings are 
used to calibrate the field dependence of the two cernox 
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sensors, while the values of ST determined with / at the 
two values provide a check on the linearity of the sam- 
ple response. By testing with a standard material that 
has no intrinsic field dependence in k xx (nylon), we have 
found that at temperatures above 8 K this method pro- 
vides an accurate and highly reproducible determination 
of the intrinsic conductance to a resolution of 1 in 10 3 . 
Although the bridge-balance method is capable of higher 
resolution, the present technique lends itself to full au- 
tomation. A higher density of points may be obtained, 
and checks (e.g. for linearity) can be made in situ. A 
pair of thermocouple junctions are used to detect the H- 
antisymmetric (Hall) gradient to obtain n xy . 

With the high sensitivity, we have found that, in un- 
twinncd, optimally-doped YBa2Cu30 x (T c = 93 K, x — 
6.95) field H || c becomes increasingly hysteretic 

below 35 K. Although the hysteresis is small (about 5% of 
the total K X x at 8 K) , it greatly complicates the extraction 
of the QP current (we discuss this later). In underdoped 
crystals, however, the hysteresis is unobservable up to 14 
tesla (less than 1CP 3 ), and the observed k xx vs. H is a 
faithful representation of its intrinsic dependence on B. 

In this report, we discuss data from a twinned, under- 
doped crystal in which T c = 63 K, and x = 6.63. The 
zero- field temperature profile of k xx is shown in Fig. [|. 
The relative magnitude of the anomaly in k xx is only 
about a quarter of that in the 93-K YBCO, but larger 
than in optimum Bi 2212 § and La 2 - x Sr x Cu0 4 (LSCO) 
0. Also shown (solid symbols) is our new estimate of 
the phonon conductivity (kb). One of our main results 
is that the entire field dependence of k xx derives from 
the QP current, while the phonon current is unaffected 
by H. 

3. Results and Analysis 

Figure displays the field dependence of k xx at selected 
T. With decreasing T, the initial slope of n xx increases 
rapidly. Below 15 K, the rapid decrease crosses over to 
an almost flat dependence, which recalls the plateau fea- 
tures observed in single-domain Bi 2212 ||. Within our 
resolution, there is no resolvable hystereses at the tem- 
peratures investigated. The higher precision and larger 
range of the new data enable us to compare the observed 
field dependence with various expressions. We find that 
the field dependence is accurately fitted (Fig. |2|) to the 
expression 



k xx (B,T) 



K e (T) 



(l+p(T)\B\ 



kb(T), 



(1) 



where the entire B dependence resides in the denom- 
inator of the first term, and the term is a field- 
independent background. At each temperature, the fit 
yields the two parameters n e (T) and p(T) associated with 
the QP current, and the term kb which we identify with 
the phonon term, viz. kb = n P h- 
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FIG. 2. Variation of k xx with field H || c in underdoped 
YBa2Cu3 06.63 at indicated temperatures. All curves are 
non-hysteretic to 1 part in 10 3 (the curves above 15 K are 
superpositions of sweep-up and sweep-down traces). A fit to 
the data at 22.5 K is also shown superposed. Curves at 15 
K and below (discrete symbols) show the approach of k xx to 
the //-independent value (kb)- 

In our previous experiments on optimally-doped YBCO 
H and LSCO @, fits to Eq. |l| were ambiguous because 
of strong hystereses (the distortions introduced caused 
the extracted p(T) to be non-monotonic in T). These 
extraneous effects led us to consider alternate expressions 
(see below), as well as a possibly field dependent n p h- 
However, the present results have clarified the problem. 
In addition to the absence of observable hysteresis, the 
curves for k xx at low T display pronounced curvature in 
moderate fields, corresponding to a strong attenuation 
of the QP current. The attenuation uncovers a nearly 
field-independent background that we identify with the 
phonon thermal conductivity. In the main panel of Fig. 
|3j, we display the T depenence of the parameters K e and 
P(T). 

The motivation for Eq. pi is that, near the vortex core, 
the steep variation of the pair potential and the circu- 
lating superfluid together present a strong scattering po- 
tential for an incident QP |[l0) . Expressing the scattering 
rate as a transport cross-section a tr , we assume additivity 
of the rates, and write the QP mean- free-path in a field 
as £(B) = 4>/[l + £o&tr\B\/4>o], where £ is the zero-field 
value of the mean-free-path, and 00 the flux quantum. 
Thus, in this model, the QP conductivity is given by Eq. 
with @,|l]| 



p(T) = £ a a tr /4>o- 



(2) 
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In the Boltzmann equation approach, we may write 
the QP thermal conductivity (in zero field) as [fl2[ 



Ke(T) 



9/ 

<9£ 



)£(k)V(k) 2 r(k), 



(3) 



where -E(k) and r(k) are, respectively, the energy and 
lifetime of a QP in the state k, and v(k) = h~ 1 \7E(k) 
its group velocity. In a <i-wave superconductor at low T, 
the excitations are confined to Dirac cones at the nodes 
where the energy may be parametrized as E(k\,k2) = 
Tiy/{kiVf) 2 + (k2V2) 2 (fci and k 2 are the components of 
k normal and parallel to the Fermi Surface), K e reduces 
to 



V k%T 2 



v 2 



-(1 + 1). 



(4) 



where £q = w/To is the mean free path at the nodes, and 
?/ = / °° dx x 3 (-df/dx) ~ 5.41. 

The T 2 dependence of K e in Eq. [| is masked by the 
strong T dependence of £q- However, in our experiment, 
the latter is obtained independently from the field de- 
pendence (with p(T) given by Eq. ||). We may di- 
vide out the T dependence of 4>> to isolate the quantity 
L e (T) = K e (T)/p(T). Comparing Eqs. | and |, we are 
left with an expression that contains only two material- 
specific parameters v 2 and at r , viz. 



, . 77 k%T 2 
7T Ti 



V 2 <T tr 



(5) 



Figure || (inset) reveals that the measured L e (T) displays 
a nearly T 2 dependence at low T, that may be fitted 
to give L e (T) = 6.9 x 1CT 3 T 2 WT/mK. Comparing this 
expression with Eq. 0, we determine from our experiment 



v 2 a tr = 2.11 x l(T 4 m 2 /s. 



(6) 



4. The Hall angle 



The field parameter p(T) has been extracted from k xx 
alone. While its nominally 1/T 2 variation (Fig. 0) is 
consistent with the identification p(T) ~ £$ (Eq. 0), it 
is important to see if this is consistent with a separate 
experiment. We turn next to the Hall conductivity n xy . 
In the previous Hall study on optimal YBCO Q], K xy was 
analyzed without the benefit of information on the diag- 
onal electronic current. From the analysis above, we may 
now extract the Hall angle tan 6qp{H) — K xy (H) / n e (H) 
as a continuous function of H at each temperature. In 
general, tan# displays strong negative curvature vs. H 
p3[ . Here, we restrict our discussion to the weak-field 
value 9qp(0). In underdoped YBCO, the small QP pop- 
ulation generates a weak thermal Hall current, and the 
uncertainties in determining 9qp are quite large (com- 
pared to optimum YBCO) [p~3|. Nevertheless, we find 
two interesting features of the Hall angle (see Fig. ||). 



First, 9qp(0) (solid triangles) and p (open circles) share 
the same T-dependence from T c to about 20 K. Secondly, 
we recall that the normal-state electrical Hall angle On 
(open triangles) follows a 1/T 2 dependence. 
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FIG. 3. The T dependence of the zero-field QP conduc- 
tivity K e and the field scale l/p(T) extracted from fits to n xx 
vs. H in underdoped YBCO. The T 2 dependence of l/p(T) 
suggests that it is related to a scattering rate. The solid line 
is co + ciT + c 2 T 2 , with co = 0.25, a =0.016, and c 2 = 0.0012 
(SI units). The inset plots the T dependence of the ratio 
L e (T) = K e (T)/p(T) (in which £0 cancels). In the vortex 
scattering model, L e should vary as T 2 in a d-wave supercon- 
ductor at low T (Eq. |) 

The new values for tan#Qp(0) lie on the curve for 9jq 
extrapolated below T c . As shown in Fig. the three 
quantities p, tan 0qp(0) and tan On fall on the same curve 
over about 2.5 decades (in the plot p and tan 6*^(0) are 
related by a constant scale factor). Just above T c , tan On 
displays a slight dip associated with fluctuation effects. 
The similarity between Oqp and On has also been pointed 
out by Zeini et al. [G3 . 

We briefly discuss our interpretation. These results 
suggest that the 1/T 2 dependence is, in fact, intrinsic 
to the QP's below T c . The ubiquitous T 2 dependence 
of cot On in the normal state appears to be an exten- 
sion of the low-temperature behavior into the normal 
state. The similarity between the T-dependences of p 
and Oqp imply that the diagonal and the Hall channels 
relax with the same T-dependence, ~ 1/T 2 , consistent 
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with simple Drude behavior. Just as in conventional 
metals, the thermal Hall resistivity in the mixed state, 
W xy = K xy /nl = tan6>Qp(0)/K e , should provide a mea- 
sure of the heat capacity of the QP population, as it does 
(see inset of Fig. ||). 
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FIG. 4. The temperature dependences of the field scale 
p(T) (open circles) and the weak-field QP Hall angle below T c 
(closed triangles) and the electrical Hall angle tanf9jv above 
(open triangles). The close similarity of T-dependences of p 
and tanf9gp strongly suggests that p is proportional to the 
QP lifetime. The broken line shows that the QP weak-field 
Hall angle is numerically equal to the extrapolation of tan On ■ 

This conventional behavior is abruptly altered when we 
cross T c into the normal state. The Hall angle continues 
to relax at the same numerical rate. By contrast, the 
scattering rate in the diagonal conductivity undergoes 
a dramatic change |l]J. The transport mean-free-path 
£q decreases abruptly by a factor of 4-6 (10 in optimal 
crystals) across T c . This sharp decrease, followed by a 
nominally T-linear scattering rate, is responsible for the 
anomalously strong T-dependence of the Hall coefficient 
in the normal state. Thus, the anomalous channel re- 
sponsible for most of the strange-metal properties is the 
diagonal conductivity. The Hall channel appears to be 
quite conventional. A detailed discussion of these results 
including optimal crystals will appear elsewhere 

5. Discussion 

We have discussed how the field dependence of 
underdoped YBCO may be analyzed to extract electronic 
parameters, such as £q and Otr- The analysis derives from 
two striking features of k xx intrinsic to high-purity 60-K 
YBCO crystals at low T, namely the steep decrease of 
k xx in weak field followed by a saturation at high field, 
and the absence of resolvable hysteresis. 



The first feature allows us to fit a much broader range 
of field scales (as expressed by the dimensionless param- 
eter p{T)B). The higher density of measurements also 
helps. We illustrate this point as follows. In a pre- 
vious attempt to analyze similar measurements in 
LSCO, it was found that the H curves were 

equally well-fitted by Eq. [l] and the expression G(B) = 
ip(l/2+B /B)+]xi(B /B), withip(x) the digamma func- 
tion. With the data scatter and the smaller range of re- 
duced fields B/Bq in LSCO, the two fits could not be dis- 
tinguished, and Ong et al. argued the case for adopting 
G(B) to describe k xx in LSCO. However, with the larger 
field scale pB here, we find that Eq. [I] provides a much 
better fit (this is evident if the two fits are compared in 
a plot versus IniJ). The physics underlying the two fit 
expressions is of course quite different. In light of the 
present work, we now favor adopting Eq. [j], instead of 
the digamma function fit, for analyzing k xx vs. H curves 
in cuprates. 

The second feature (no observable hysteresis) relates to 
the issue of remanence and vortex pinning in cuprates. It 
is known that the relaxation of non-equilibrium flux dis- 
tributions may produce a slow drift in f^xx 

if it is observed 

a few seconds after a change in H is made ||||. Further, 
vortex pinning effects at low T can lead to step-like jumps 
when the field sweep direction is changed. In opti- 
mum YBCO and in strongly overdoped Bi 2212, we find 
that k xx increases step-wise when the sweep direction of 
H is reversed from up to down. Recently, however, a 
hysteretic loop of the opposite sign has been reported by 
Aubin et al. || in Bi 2212 (k xx decreases step-wise when 
H is swept down). At present, the origin of the hystereses 
in k xx is not understood (especially the existence of hys- 
tereses with different signs). We note that the magnitude 
of the hystereses in is much smaller than that observed 
in the magnetization M vs. H. In our measurements on 
underdoped YBCO, no hysteresis in f^xx 

observed for 

fields as large as 14 T at temperatures down to 6 K, even 
though hystereses are sizeable in the M vs. H curves. In 
particular, the magnitude of k xx at the plateau-like re- 
gion is not hysteretic. The absence of hysteresis implies 
that the magnetization is too small to influence the mea- 
sured k xx , so we may assume that B — fioH, as tacitly 
assumed in the fits. 

By contrast, hysteretic effects cannot be neglected in 
optimally doped YBCO (as discussed above). Stronger 
vortex pinning is clearly responsible for the larger hys- 
teresis in the 93-K crystal. Below 35 K in this crystal, 
the hysteresis steadily increases to about 5% at 8 K. Al- 
though the hysteresis is small, the remanence produces 
in the trace of 

f^xx VS. H both a broadening at small 
H and an asymmetry about H = that strongly distort 
the fit to Eq. [l| The distortions preclude a meaning- 
ful extraction of below about 35 K. Thus, of the various 
phases of the cuprates we investigated (YBCO, Bi 2212 
and LSCO), the underdoped phase of YBCO appears to 
be the most suitable for our purpose of isolating the QP 
current from the total thermal current. 
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In high-purity single-domain crystals of Bi 2212, the 
field dependence of k xx displays a distinct break in slope 
in k xx at a characteristic field Hk , followed by a plateau 
region in which n xx is nearly independent of H |5|. The 
field Hk, which varies approximately as T 2 , was inter- 
preted as a field-induced phase transition, possibly in- 
volving a new order parameter. We compare the present 
results with the two findings in Bi 2212, i.e. the kink fea- 
ture at Hk and the existence of the plateau. As shown in 
Fig. ||, k xx in 60-K YBCO smoothly crosses over into the 
field-independent region at low T, instead of displaying 
a sharp kink. While the plateau regime is similar in the 
two systems, the kink feature signalling a phase transi- 
tion is absent in the YBCO crystal. A difference between 
the two systems is the electronic anisotropy. From the 
resistivity anisotropy p c /p a b (~ 10 5 compared with 10 3 ), 
Bi 2212 is much closer to the 2D limit than underdoped 
YBCO. Whether this is a significant factor is a subject 
for future investigation. 

In their experiment Aubin et al. (§] observed the value 
of k xx in Bi 2212 at the plateau (at 8 K) to be 1% higher 
in the field sweep-up direction than in sweep-down. In 
rough analogy with the magnetization profile in the Bean 
model, they raise the issue that the plateau may be as- 
sociated with, or reflect a specific state of the vortex sys- 
tem. In our response ||, we pointed out that the hys- 
teresis in their sample is about 5 times larger (at 8 K) 
than in the two crystals used by Krishana et al. ||. At 
higher temperatures, 15-20 K, where the plateau is just 
as prominent, the hysteresis is almost unresolved. The 
hysteresis is an extrinsic effect possibly associated with 
stronger flux pinning in a more disordered crystal. The 
present results lend a fresh perspective to the question 
whether the plateau is intrinsic to the QP system or the 
product of a particular state of the vortex system. The 
YBCO results show that, whenever the QP current is 
very strongly suppressed by the available field, the ther- 
mal conductivity that remains is indeed field-independent 
and non-hysteretic. This shows that, in YBCO, a plateau 
regime definitely exists; the lack of observable hysteresis 
shows that it has nothing to do with a particular state of 
the vortex system. However, to access it, it may be nec- 
essary to work in the underdoped regime and to use high- 
purity crystals with weak pinning (as discussed above, we 
are unable to access the plateau region in 90-K YBCO). 

In principle, the analysis described yields quantitative 
information on the quasi-particles. Having both the di- 
agonal and off-diagonal conductivities available reduces 
the uncertainties in identifying the measured parame- 
ters, as well as provides consistency checks. The weak- 
field Hall angle may be expressed as a 'Hall' mean-free- 
path £h = 0QphkF /e. The value of 8qp at 10 K gives 
Ih ~ 4, 200A The similarity of the T dependences in 
6qp and p implies that Ih is proportional to £q. If the 
proportionality constant is 1, we may use In in Eq. |^ to 
find that <7tr ~ 90A, about 1.7 times the diameter of the 
vortex core (2£) (£ ~ 26 A if the upper critical field H C 2 ~ 



50 T). Using this value of a tr in Eq. |L we obtain the ve- 
locity v 2 ~ 2.3 x 10 6 cm/s. From the penetration depth 
variation AA = 4.3 A/K Q, Wen and Lee 0] obtain 
the velocity anisotropy Vfjvi ~ 7.6. With our estimate 
for V2, we find that the Fermi velocity Vf ~ 1.8 x 10 7 
cm/s. 
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